Tunneling-driven breakdown of the 331 state and the emergent Pfaffian and 

composite Fermi liquid phases 
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We examine the possibility of creating the Moore- Read Pfaffian in the lowest Landau level when 
the multicomponent Halperin 331 state (believed to describe quantum Hall bilayers and wide quan- 
tum wells at the filling factor v = 1/2) is destroyed by the increase of tunneling. Using exact 
diagonalization of the bilayer Hamiltonian with short-range and long-range (Coulomb) interactions 
in spherical and periodic rectangular geometries, we establish that tunneling is a perturbation that 
drives the 331 state into a compressible composite Fermi liquid, with the possibility for an inter- 
mediate critical state that is reminiscent of the Moore-Read Pfaffian. These results are interpreted 
in the two-component BCS model for Cauchy pairing with a tunneling constraint. We comment on 
the conditions to be imposed on a system with fluctuating density in order to achieve the stable 
Pfaffian phase. 

PACS numbers: 73.43.Cd, 73.21. Fg, 71.10.Pm 



I. INTRODUCTION 

When electrons are confined to a two-dimensional (2D) 
plane and subject to a strong perpendicular magnetic 
field, they organize themselves into fascinating strongly- 
correlated quantum phases. The most prominent exam- 
ples are the Laughlin states^ = $>2k+i({z}), which 
may be written in terms of the Laughlin- Jastrow factor 
$ m ({z}) = l\i<j( z i - z j) m , where Zj = Xj + iyj denotes 
the complex coordinate of the j-th electron, m, k are in- 
tegers and we have neglected the Gaussian factor ubiq- 
uitous in the lowest Landau level (LLL). Laughlin states 
describe the fractional quantum Hall effect (FQHE) that 
occurs when the filling factor v = N/N^, which is de- 
fined as the ratio between the number of electrons N 
and the number of flux quanta Nj, threading the 2D 
system, is a simple fraction with an odd denominator, 
v = l/(2k + 1)3 Laughlin's construction has been gen- 
eralized, within the framework of the composite-fcrmion 
(CF) theory— to explain the rich phenomenology of a 
whole sequence of observed odd denominator fractions. 
According to the CF theory, one understands the FQHE 
as an integer quantum Hall effect in an effective mag- 
netic field that vanishes at v = 1/2. As a consequence, 
CFs then form a compressible CF Fermi liquid (CFL)^£ 
|*fl) = Vlll dct [4n(zj)]<&2{{z}), as seen in the absence 
of the Hall plateau in single thin layers at v = 1/2 in the 
LLL<£ Here, as before, the characteristic quantum Hall 
correlations are captured in the Jastrow factor "^({z}) 
which we refer to as the charge part of the wave function 
(as it carries the flux through the system) and 4>i{zj) 
are the single particle states (an overall projection Vlll 
to the LLL may be needed to yield a proper trial wave 
function that is analytic in Zj.). 

However, some quantized even denominator states 



exislii and are usually associated with the first excited 
Landau level, where the nature of the effective interac- 
tion is believed to facilitate the pairing between CFs>£ 
The paradigm of such paired states is the so-called "Pfaf- 
fian" state introduced by Moore and Read,— |\T/pf) = 
Pf [l/(zi- Zj)]$ 2 ({z}), which explains the FQHE ob- 
served at v = 5/2 3 In addition to the charge part 
<3>2({z}), which fixes the filling factor, we have also a 
pairing in the neutral sector described by the object Pf M> 
In contrast to the Laughlin and Jain states with any- 
onic excitations satisfying Abelian statistics, the Moore- 
Read state represents the simplest paired state of spin- 
polarized electrons which supports excitations with non- 
Abelian statistics^ of interest in topologically protected 
quantum computation^ 

If the spin of the electrons is not necessarily frozen 
out by the magnetic field, the electrons could find it 
more favorable to reorganize themselves into one of 
the competing Abelian phases called multicomponent 
Halperin states*^ In these states the Hall quantiza- 
tion is a result of internal degrees of freedom of the 
electrons (provided by the spin or layer index). At 
half filling, a two-component candidate is the 331 state, 

|*33i> = *3* ra ({«t})*g ,trB ({^})*? ,ter ({-»r,«4.}). writtcn 

as a straightforward generalization of the Laughlin state 
to two species of electrons f an d I- In order to sat- 
isfy the constraint of the fixed filling factor v = 1/2, 
apart from the usual Laughlin- Jastrow factors between 
the electrons of the same species ($ intra ), one also has 
to account for the inter-species correlations through the 
factor $)f ter ({z, w}) = ILj^ — w j) k - Alternatively, us- 
ing the Cauchy identityi 3 - we can cast the 331 state into 
the form |* 33 i) = det [l/(z,:,t - Zj t ±)} § 2 {{Z}) which ex- 
tracts the charge part $2 (Z's denote all particles regard- 
less of their spin index) from the neutral part where the 
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pairing is described in terms of a Cauchy determinant 
between f and | particles. Numerical calculations 1 ^— 
indicate that Halpcrin's 331 wave function is likely to be 
at the origin of the v — 1/2 FQHE in bilayer quantum 
Hall system a 19 ' 20 as well as in wide quantum wells £L 

In this paper we investigate whether it is possible to 
create the Moore-Read Pfaffian in the LLL by convert- 
ing the two-component 331 state into a single-component 
state. Mathematically, this is easily achieved by antisym- 
metrizing the neutral (Cauchy determinant) part of the 
1^331) between j" and j. . 22 i 23 However, such a procedure 
is a very complex mathematical entity because it creates 
a state with different physical properties (non-Abelian 
statistics out of the Abelian), while we are interested in 
a physical mechanism that mimics the antisymmetriza- 
tion in an experimental situation. We restrict the discus- 
sion to the Coulomb bilayer system, which is a generic 
two-component system where the "spin" a =7, |. denotes 
the two layers in which the electrons are localized. In 
such a system, it is commonly speculated that the an- 
tisymmetrization mechanism is provided by the tunnel- 
ing term ~ — AsasS x , which combines the single-particle 
wave functions into symmetric (even) and antisymmetric 
(odd) superpositions, | ± |. The tunneling term favors 
the even superposition (channel) where one expects to 
find a weakly paired (Moore-Read) phase. We establish 
that this route towards the Moore-Read state is compli- 
cated by the presence of the compressible CFL, which is 
the resulting phase for large tunneling. Along the way, 
one may arrive at a critical state that shares some prop- 
erties with the Moore-Read Pfaffian, but we do not find 
evidence that this state represents a stable phase. Re- 
cent experiment o 24 ' 25 found even denominator fractions 
in wide quantum wells in the LLL. The results obtained 
within the present bilayer model may be relevant also in 
the study of wide quantum wells insofar as the latter can, 
with moderate approximations, effectively be described 
by a bilayer Hamiltonian, where the tunneling term mim- 
ics the effective confinement gap between the lowest and 
the first excited electronic subleveli^Z 

The remainder of this paper is organized as follows. In 
Sec. [II] we introduce the BCS model for spinful fcrmions 
with tunneling, first proposed by Read and Green^ This 
model decouples into an even and odd channel, and the 
tunneling term increases the population of the even chan- 
nel, where one expects to find the weak pairing (Pfaffian) 
phase. However, the increase of tunneling also leads to 
the effective weakening of the coupling, which eventu- 
ally drives the system into a CFL phase when the total 
number of particles is held fixed. This is the situation we 
analyze in detail in our exact diagonalization calculations 
in spherical and toroidal geometries (Sec. IHIj) . These re- 
sults are discussed in the context of the phase diagram 
of Ref. (see Sec. MI C[) . We furthermore introduce a 
generalized tunneling constraint which, in the BCS de- 
scription, leads to a stable weak pairing (Pfaffian) phase 
in the even channel when the density of the of the sys- 
tem is not fixed (Sec. IIV|) . We present our conclusions 



in SecEl 



II. BCS MODEL WITH TUNNELING 

At v = 1/2, the CFs experience a zero net magnetic 
fields and if we limit ourselves to the neutral part of 
I "Jpf ) , they may be described within the framework of the 
effective BCS model introduced in Ref. H^. We consider 
the system to be at zero temperature and neglect fluctu- 
ations in the Chern-Simons gauge field that are related 
to the charge part of l^pf).— The Hamiltonian which 
describes the Cauchy pairing between 'f and 4- particles 
with tunneling As^s reads 

H = E [ e ~k(4t c kt + 4i c n) + (Ak4 t ci H + h.c.) 

k 

-A^(4 t c kl + C ki c kf) ^ ^ 

where ?k = ek — fx, in terms of the putative CF disper- 
sion relation ek and the chemical potential fi, which is 
assumed positive fi > 0. Notice that because of the van- 
ishing net magnetic field, the 2D wave vector k = (k x , k y ) 
is again a good quantum number. The order parameter 
Ak = Ao(k x — iky) is chosen to describe p-wave pairing, 
and we assume that Ak and [i are not renormalized by 
the tunneling. 

With the help of the even, Ck, e = (ck,t + c k,j.) /V2, and 
odd spin combinations Ck, D = (ck,t — C k.,i)/V2~, the Hamil- 
tonian (fTJ) decouples into an even and odd channel^ 
H = H e + H°, where (the index r denotes the even and 
odd channel, r = e, 0) 

RT = E [( £ k - f)cl T C K r + (A k C kiT c! kir + H.C.)] , 
k 

(2) 

in terms of the chemical potentials /i e = /i + Asas/2 
and n° = fx — As as/ 2 for the even and odd channels, 
respectively. Furthermore, the even/odd p-wave order 
parameters read A k = Ak/2 = (Ao/2)(k x — ik y ) and 
A° = A_ k /2 = -(Ao/2)(fc a - iky). 

For moderate tunnelings, the effective chemical po- 
tential /i cff of the whole system may be viewed as the 
weighted sum of the two channels, /i cff = Pfi e + (1 — P)(J,°, 
where P measures the population of the even channel 
(1/2 < P < 1) and may have a complicated dependence 
on Asyjg. In particular, for some values of As as we 
may be below the critical line /J, cff = PAsas and inside 
the non-Abelian (Pfaffian) phase. However, in the limit 
of large tunneling, the system is dominated by the even 
channel and the chemical potential of the whole system is 
because P = 1 . Remember that the associated 
BCS wave function in the even channel reads 

\iPbcs) = ]J i 1 + f kC ik,e c -k :e ) Ivacuum), (3) 

k 
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in terms of the pairing function ~ Wk/""k ~ A* e /Ao. 
One notices then that an increase of the chemical po- 
tential /i e controlled by the large value of the tunneling 
parameter As^s is equivalent to a reduction of the order 
parameter Ao. Therefore the BCS system will eventu- 
ally be transformed into the one of the Fermi liquid. We 
can see this more explicitly by examining the relevant 
excitations of the even channel^ 



also consider long-range Coulomb interaction, 
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in the limit of large [f around k = |k| = 0. They be- 
come unstable and k = becomes a point of local max- 
imum. The minimum is expected to move to |k| = kp, 
the Fermi momentum.— Therefore if Ao does not "renor- 
malize" with increasing As^s, the net effect of the strong 
tunneling (fi e ^> fi) on the Cauchy pairing is to drive the 
system into a Fermi liquid. This is not unexpected be- 
cause one retrieves an effective one-component system in 
this limit, where all particles are "polarized" in the even 
channel. The Pfafhan physics may however play a role 
in the intermediate state before complete polarization. 
We revisit the BCS approach in Sec. IIV1 with a slightly 
different perspective in which the antisymmetrization is 
imposed, in a functional formalism, with the help of a 
Lagrangian multiplier which plays a similar role as the 
present tunneling term As^g. 

As we pointed out earlier, the population of the even 
channel P may be a complicated function of tunneling. 
In the following section we use exact diagonalization of 
small finite systems in order to get a hint of the form of 
this dependence P = P(Asas) and determine the nature 
of possible phases as P increases from 1/2 to 1. 



III. EXACT DIAGONALIZATION 

Here we study the full interacting quantum Hall bilayer 
Hamiltonian for small finite systems in the presence of 
tunneling ; 16 ' 18 
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where in coordinate representation we have 25*3; = 

/ dr \I'|(r)^'4_(r) + H.c. , ^(r) creates a particle at the 

position r in the layer a. We have decomposed the in- 
teraction into terms between electrons belonging to the 
same layer (y lntra ) and those residing in opposite layers 
(V lnter ). We consider a short-range interaction, defined 
as 

VH?*{r) = FxV 2 5(r), V£?*(r) = V 5(r), (6) 

which produces the 331 state as the densest and unique 
zero energy state when Vq, V\ are chosen positive! 10 ' 28 We 



^oufM = e 2 /er, 



'Colli 



(r) = e 2 /evV+d 2 , (7) 



where d is the distance between layers. We fix the to- 
tal number of particles in our calculations to be an even 
integer and, unless stated otherwise, take d = Ib (Ib is 
the magnetic length) , which merely sets the range for the 
distance between the layers where the Coulomb ground 
state is supposed to be fairly well described by the 331 
wave function. Confining the electrons to a compact sur- 
face such as the sphere^ 8 - or torus^ the Hilbert space 
becomes finite and one may exactly diagonalize the inter- 
acting Hamiltonian ([5]).— The ground state obtained in 
this way can be numerically compared with the trial wave 
functions 1^331) and |\&pf) as a simple scalar product be- 
tween vectors in the Hilbert space. In these calculations 
I^Pf) is defined in the even basis i.e. single-particle states 
are understood to be even combinations of the original 
t, I states. 



A. Sphere 

If we wrap the electron sheet into a sphere and place 
a magnetic monopole in the center which generates 
magnetic flux quanta perpendicular to the surface, we are 
left with a finite basis of single particle states indexed by 
0, N$. Translational symmetry in the plane is replaced 
by rotational symmetry, which leads to a classification 
of the many-body states by the eigenvalues of angular 
momentum L and its z-component L z .— The two-body 
interaction such as (JSJ) and (O is parametrized by Hal- 
dane pseudopotentials V[ a which represent the energy 
of a pair of particles located in layers a, a' with relative 
angular momentum Incompressible quantum Hall 
states are invariably obtained in the L = sector of the 
Hilbert space. They are further characterized by a topo- 
logical number called the shifts 5, defined through the 
relation = — 5. To specify uniquely the Hilbert 
space corresponding to the given trial state, one needs to 
know S to get the correct value for the pair of (iV, A^,). 
Sometimes it is possible to find different values of v, 5 
that yield the same (TV, N$) - this is called the aliasing 
problem because two different trial states get realized in 
the same Hilbert space. We disregard such cases in our 
calculations. 

In FigfT]we present results of exact diagonalization on 
the sphere for the short-range (331) Hamiltonian (j6]) and 
long-range (Coulomb) Hamiltonian ([7]). 331 and Pfaf- 
fian trial states occur at the same value of the shift, 
thus we are able to track their evolution as a function 
of tunneling simultaneously. We also use (S x ), the ex- 
pectation value of the S x operator in the ground state, 
to monitor the two-component to one-component transi- 
tion, whereas (S z ) remains zero throughout, which is due 
to the weaker interlayer as compared to the intralayer 
interaction. Starting from the long-range Hamiltonian 
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(Fig. [T] left panel), we see that the 331 state gives way 
to a Pfaffian-like ground state, with the overlap quickly 
saturating to a value of around 0.92 . 16 ' 18 The transition 
occurs for Ag^g ~ 0.1e 2 /ds which agrees well with the 
typical experimental value and shows little size depen- 
dence when the largest accessible system N = 10 is con- 
sidered (note that the subsequent N — 12 system suf- 
fers from the aliasing problem). On the other hand, no- 
tice that for the short-range Hamiltonian (Fig. [T] right 
panel), the 331 state is much more robust to perturba- 
tion by A5^5: before it reaches full polarization in the in- 
direction (maximum (S x )), the overlap with both incom- 
pressible states drops precipitously beyond some critical 
Asas which is rather size-sensitive (it also depends on 
the values of the parameters one chooses in Eq. ([BJ , but 
the qualitative features of the transition are reproduced 
for many different choices of Vq , V\ ) . 
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FIG. 1: (Color online) Overlaps between the exact ground 
state of the Coulomb bilayer (left panel) and short range 331 
Hamiltonian (right panel) with the 331 state (O331) and the 
Pfaffian (Opt), as a function of tunneling Asas- Also shown 
on the right axis is the expectation value of the S x component 
of pseudospin (for N = 8 system) which characterizes the two- 
component to one-component transition. 

We see nonetheless that the breakdown of a two- 
component phase yields a one-component state mani- 
fested by (S x ) -> N/2 (this is the limit P ->■ 1 from 
Sec|TT|. At the transition, (S x ) develops a small but vis- 
ible kink. Focusing on the large tunneling limit, we find 
that the nature of the ground state is effectively that of 
the single layer (polarized) ground state for the symmet- 
ric interaction V+(r) = [V intra (r) + V intcl '(r)] /2. This 
intuitive result was directly verified for all the available 
system sizes, including very large A = 10 system on both 
sphere and torus (see Sec. IIII B[) . In view of this, it is 
not surprising that the large-tunneling limit of the short- 
range Hamiltonian is the compressible CFL: V + in this 
case reduces to the repulsive hard core V\ pseudopoten- 
tial which has a tendency to produce the Jain CF state. 
This is also apparent in the fact that the ground states for 
large tunneling are obtained in the angular momentum 
sectors that agree with the predictions for the excitations 
of the CFL yielding overlap of 0.99 with the excited CF 



sea ground state. 

Therefore, the results for the short-range Hamiltonian 
are suggestive that wc may have a direct 331-CFL tran- 
sition in the thermodynamic limit because the transition 
point seems to be shifting towards smaller tunnelings as 
we increase N . Notice, however, that in contrast to the 
incompressible 331 and Pfaffian states, which occur at a 
shift 5 = 3, the CFL has a shift S = 4. On the sphere, 
the two incompressible states can thus not be directly 
compared to the CFL, and the evidence for the 331-CFL 
transition is therefore indirect. This problem is circum- 
vented in ED on the torus presented hereafter in Sec. 
IIII Bl In the Coulomb case, on the other hand, we ob- 
serve a curious saturation of the ground state overlap 
with the Pfaffian. Wc attribute this feature to the effect 
of the long range Coulomb potential on a finite system. 
One notices that by adding an asymptotic tail to the 
"intra" component of the short-range pscudopotentials 
ymtm = ymtm + a / 2y /Z (a nonzero for L > 3) , one 
progressively increases the critical value of Asas for the 
abrupt drop of the overlaps as a — > 1 (pure Coulomb) . In 
fact, for A~ = 8, 10 it is sufficient to consider only V^ ntra 
to achieve the saturation and push the critical value of 
A SA s to infinity. 

Results for the Coulomb interaction in the large Asa.i 
limit (Fig. [IJ are similar to those obtained in Refl3ll 
where single-layer Zhang-Das Sarma interaction was 
used. As long as we are in the large Ag^g limit, V + (r) 
interaction produces numerically the same effect as the 
Zhang-Das Sarma interaction. In particular, transition 
to a Moore-Read Pfaffian will be induced if the layer 
separation d is sufficiently large^ Of course, these two 
interactions arc different from each other and the fact 
that they yield the same phenomenology (phase transi- 
tions as d is varied) only means we are probing a critical 
state where even the slightest perturbation away from 
pure Coulomb interaction (coupled with the bias of the 
shift) is sufficient to generate incompressibility. However, 
despite large overlap, the gap remains very small after the 
transition. The difference between the two interactions is 
obvious in the torus geometry (Scc lHI B"j) . The new result 
of the present paper is that we find the Pfaffian signature 
even in the region without full polarization (P < 1), as 
we elaborate in Sec lIII CI 



B. Torus 

Another way to compactify an infinite plane is to im- 
pose periodic boundary conditions on a unit cell a x 6i 28 ' 29 
This also produces a finite set of N$ = ab/2nl 2 B single- 
particle states which are periodic functions under the 
transformation of coordinates x — > x + a,y — > y + b 
(we assume that there are no additional phases gener- 
ated by such a transformation). Because of the pres- 
ence of magnetic field, the many-body Hamiltonian re- 
duces in the invariant subspace of the magnetic trans- 
lation group^ and the eigenstates are labeled by the 
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pscudomomcntum k, restricted to the Brillouin zone 
(2-Ks/a, 2irt/b); s, t = 0,...,N — 1, where N is the great- 
est common divisor of N,N^. Contrary to the sphere, 
trial states on the torus are uniquely specified by their 
filling factor and thus we can directly address transitions 
between the 331 state, the Pfaffian and the CFL. More- 
over, states in this geometry are distinguished by their 
ground-state degeneracy. If the filling factor is v = p/q 
(p, q having no common divisor) , there is a generic degen- 
eracy of q which comes from the center of mass motion^ 
having no physical importance, so we implicitly assume 
it in what follows. Apart from this, there can be ad- 
ditional degeneracies occuring due to special point sym- 
metry of the Brillouin zone (trivial) or those that arise 
from the multicomponenlj^ or the non-Abelian nature of 
the statei&i^ For the 331 state we expect a quadruplet of 
states (up to the center of mass degeneracy) one of which 
belongs to the k = sector of the Hilbert space and the 
remaining three are located at the corners of the Brillouin 
zone, k = (0,N/2), (N/2,0), (N/2, N/2). In contrast to 
the 331 state, the Moore-Read Pfaffian has only a three- 
fold degeneracy^ k = (0, N/2), (N/2, 0), (N/2, N/2), 
whereas compressible states in general do not possess 
clearly defined degeneracies (they may appear to have 
accidental degeneracies which are functions of the aspect 
ratio of the torus a/b, particle number and any other 
parameter). These are the expectations based on the 
analytic form of the trial wave functions and their par- 
ent conformal field theories but they are also borne 
out exactly in the numerical diagonalization of the model 
Hamiltoniansi&2i 
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FIG. 2: (Color online) Low-energy part of the spectrum (rel- 
ative to the ground state) of the short-range 331 (left panel) 
and the long-range Coulomb Hamiltonian (right panel) for 
TV = 8 electrons on the torus at v — 1/2 and aspect ratio 
a/b — 0.97 as a function of tunneling Asas (right axis). Also 
shown (left axis) is the overlap with the phases we identify as 
the 331 state and the CFL. 

In Fig|2] we show the relevant low energy part of the 
spectrum of the Hamiltonians and (JT]) as a function of 
tunneling, measured relative to the ground state (right 
axis), for TV = 8 electrons and the fixed aspect ratio 



0.97 in the vicinity of the square unit cell. We identify 
the multiplet of 4 states that build up the 331 phase, 
whose exact degeneracy for the short-range Hamiltonian 
and small As as (left panel) is partially lifted in case of 
the Coulomb interaction (right panel). 331 phase is de- 
stroyed for sufficiently large As as when the k = (1,1) 
state (fourfold degenerate) comes down and eventually 
forms a gapless branch with (0,4) and (4,0) members 
of the 331 multiplet (other excited states not shown). 
We identify the large-tunneling phase as the CFL phase 
because exactly the same spectrum is seen in a single 
layer with Coulomb interaction and the same aspect ra- 
tio. This transition is quantitatively reflected also in the 
overlap with the trial 331 states and CFL as a function 
of tunneling (Fig. [51 left axis). In order to take into ac- 
count the ground state degeneracy, we define overlap on 
the torus in the following way: 

|<*trial|*>| 2 = n^-T £ |(* trial (k)|vl/(k))| 2 , (8) 



trial 



keStria 



where 6>t r i a i stands for the degenerate subspace expected 
for | Atrial) • This amounts to adding together the overlap 
squared for each of the expected members in the ground 
state multiplet (normalizing the sum by the expected 
ground state degeneracy |<5>triai| to be 1 at maximum) 
and the definition is obviously meaningful only in the 
case where we had previously established the correct level 
ordering in the spectrum. 

Upon a closer look at Fig. [5J one notices that the torus 
spectra suggest little qualitative difference between the 
short-range and the Coulomb Hamiltonian. In particu- 
lar, we do not see any indication of the Pfaffian threefold 
ground-state degeneracy for large Asas which could be 
expected from the large overlap on the sphere (Fig. [TJ. 
To reconcile these two results, we again focus on the large 
tunneling limit and vary the aspect ratio of the torus to 
investigate the possibility of an emergent Pfaffian phase 
(Fig. [3]). We assume that in the large tunneling limit, we 
have effectively a single layer (polarized) ground state for 
the symmetric interaction V + (r). In Fig. [3] we show the 
spectrum of the single layer system of IV = 14 electrons 
interacting with V + (r) as a function of aspect ratio and 
connect the levels that have the quantum numbers of the 
Moore-Read Pfaffian. We also include the background 
charge correction^ One notices that, with the exception 
of a very narrow range of aspect ratios around 0.4, there 
is no evidence of a clear Moore-Read degeneracy. A nar- 
row region where we see the threefold multiplet of states 
for N — 14 also exists for N = 8, but is obscured by the 
presence of higher energy levels in systems of N = 10 and 
12 electrons. Thus we conclude that it cannot represent 
a stable phase, but a possibility remains that it is a criti- 
cal phase which becomes stronger as one approaches the 
thermodynamic limit or as one changes the interaction 
away from the pure Coulomb. 

We note that varying d (at the fixed aspect ratio) does 
not lead to any qualitative change in the ground state 
degeneracy as long as V + (r) interaction is used. This is 
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Aspect ratio 

FIG. 3: (Color online) Spectrum of the N — 14 electrons in 
a single layer at v = 1/2 interacting with V + (r) (d = 
as a function of aspect ratio. We highlight the states with 
quantum numbers of the Moore-Read Pfaffian. 

clearly different from Zhang-Das Sarma interaction which 
induces level crossings in the spectrum in such a way that 
for large d (typically beyond 41 g), a Pfaffian degeneracy 
is seen for big enough systems^ with the exception of 
N = 10. This is similar to the results in the second Lan- 
dau level^, as well as the calculations on the sphere^, 
but the prohibitively small gap suggests that such a state, 
if it exists, is very fragile. 

C. Pfaffian signatures for intermediate tunneling 
and a proposal for the phase diagram 

We conclude this section with a summary of our exact- 
diagonalization results in the two geometries in order to 
make a connection with the BCS analysis of Sec. [XT] and 
sketch possible paths of the v = 1/2 twocomponent sys- 
tem with tunneling in the phase diagram of Read and 
Green^ see Fig. |U In Fig. 2J [i has the meaning of 
the effective chemical potential /i off of the whole system 
as in Sec. [Til renormalized by Asas, he. /i = /i(Asas)- 
It is assumed that it can be approximated by the value 
of the chemical potential of the dominant even channel, 
/i cff ~ /i e and the separation between the Abelian and 
non-Abelian phases in Fig. 2] is defined by setting then 
the value of the chemical potential of the odd channel 
to zero, i.e. fj,° ~ ^i(Asas) — Asas = 0. This approxi- 
mation renders necessary taking into account the renor- 
malization of the parameters in the BCS Hamiltonian ([T]) 
with tunneling, as in Ref. [26l 

On the sphere, we first recall a very large difference 
in A saS' t ne critical value of tunneling required to fully 
polarize the system in the x-direction, for the two interac- 
tions considered. A much larger value for the short-range 
331 Hamiltonian suggests that the chemical potential for 
the even channel in this case is much more strongly renor- 
malized than for the long-range Coulomb interaction and 
therefore such a system may assume a phase trajectory 




FIG. 4: Possible outcomes of tunneling on a twocomponent 
system such as the transition to a Fermi liquid (a) or to a 
critical Moore-Read Pfaffian (61,62), in the context of phase 
diagram after Read and Green.— Note that the value for fj,p ~ 
1/m* is interaction-dependent due to the renormalized CF 
mass m* and we may have different dividing lines fj, = /j,f 
depending on the kind of the interaction.— 



labeled (a) in Fig. 0] directly moving from 331 state 
through the Abelian phase and into a CFL. 

A question we ask at this point is whether the CFL, 
a likely phase at P = 1, leaves room for other one- 
component states to form as we increase the tunneling. 
In particular, is there a possibility for a system to evolve 
along the trajectory which terminates at (61) or touches 
(62) the critical line that separates the Abelian from the 
non-Abelian phase in Fig. |4p Such an intermediate phase 
could possess significant overlap with the Moore-Read 
Pfaffian, but it would necessarily have a small gap and 
we refer to it as "critical Pfaffian" . 

On the sphere, a suitable system to detect the sig- 
nature of the critical Pfaffian is the Coulomb N = 10 
system where the large-tunneling phase is compressible 
for d < 0.51b ^ We therefore fix d — OAIb and vary 
as (Fig. O. For As as = 0, we are still largely in 
the 331 phase and for large Asas we are m the CFL; 
however, for intermediate tunnelings we see a develop- 
ing Pfaffian that establishes in a narrow range around 
Asas = 0.04e 2 /e/s- Therefore, despite "weaker" incom- 
pressibility for small Ag^s and full compressibility for 
large Ag^g, for intermediate tunneling we find evidence 
for the Pfaffian, as suggested by the trajectory 62 in Fig. 
3J. For larger values of d, as we remarked in Sec. IIII Al 
we observe saturation of the overlap for Coulomb inter- 
action and the tentative trajectory in that case resembles 
61. 

The effects of CFL physics are rendered more trans- 
parent in the torus geometry, where we have identified 
the dominant phases as 331 and CFL (Fig. [5]), with a 
direct transition between the two of them. We choose 
a value for Asas = 0.03e 2 /e^ which places the system 
in the center of the transition region (compare also with 
Fig. [5]) and examine the spectrum of an ./V = 8 system 
as a function of aspect ratio for an emerging Pfaffian de- 
generacy, Fig. [5J In torus geometry there is no subtle 
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FIG. 5: (Color online) Overlaps between the exact ground 
state of N — 10 electrons on the sphere with the 331 state 
(O331) and the Pfaffian (Opt), as a function of tunneling 
Asas, for the Coulomb bilayer Hamiltonian and d = QAIb- 

dependence on d, so we take as before d — Ib- In agree- 
ment with the results on sphere, we find a region of aspect 
ratios where the correct Pfaffian degeneracy is visible. 
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FIG. 6: (Color online) Low-energy part of the spectrum (rel- 
ative to the ground state) of the Coulomb Hamiltonian (with 
d — Ib) for N — 8 electrons on the torus at v = 1/2 and 
Asas = 0.03e 2 /e/s, as a function of the aspect ratio. Shaded 
region represents the tentative phase with the Pfaffian degen- 
eracy. 

Previous results lend support to the scenario of an in- 
termediate critical phase in a long-range Coulomb sys- 
tem, which has a small gap (Fig. [(J) but possesses large 
overlap with the Moore- Read Pfaffian (Fig. [5J. We 
would like to stress that all of our conclusions are based 
on the idealized bilayer Hamiltonian (|SJ). As such, it is 
not clear at present to what extent they apply to the 
experiment a 24 ' 25 where the electron density-profile dif- 
fers significantly from an ideal bilayer. With respect to 
theoretical studies, stronger indication of topological de- 
generacy is likely to be found in a model that assumes 
non-zero thickness of each layer^ but that would lead 
also to a substantial decrease of the gap<2£ 



IV. GENERALIZED TUNNELING 
CONSTRAINT 

In Sec. IIIII we found that in a system with a fixed 
number of particles and the tunneling commonly ex- 
pressed as — AsasSx, there is only circumstantial evi- 
dence for the clear Pfaffian phase in finite systems that 
can be studied by ED. This evidence appears most strik- 
ing when Coulomb overlaps in the spherical geometry are 
considered (Fig. [I). However, these overlaps must be 
interpreted with caution: on the sphere we can indeed 
only study the competition between Pfaffian and an ex- 
cited CFL state (containing a quasiparticle) because the 
ground-state candidates occur at different values for the 
shift 5 as mentioned in Sec. IIII A[ As we have explained 
there, the latter state is indeed favored by short-range 
interactions. Moreover, even when the ground state is 
described by the Pfaffian state, the energy spectrum is 
not that of a typical incompressible state with a gap to 
all excitations. Therefore, an explanation for the large 
Coulomb overlap with the Pfaffian may be the shift which 
artifically favors it over the CFL phase. This interpreta- 
tion agrees with the results in the torus geometry, which 
treats both phases on the same footing and which sug- 
gests that the CFL is the likely outcome of tunneling on 
the 331 state. We can give two general arguments for this. 
First, strong tunneling has a tendency to polarize the spin 
in the x-direction, and one therefore crosses over to an ef- 
fective one-component system that in the LLL favors the 
formation of a compressible CFL phase. Second, we have 
shown in the BCS approach for the charge-neutral sector 
(Sec. |H]) that tunneling docs not only lower the chemical 
potential fj,° in the odd channel but also increases that in 
the even channel, [f . We argued that this leads to the 
insufficiency of the BCS model description, which then 
describes a local maximum and the system crosses over 
to a CFL. 

It is clear that, in addition to tunneling, one also must 
find a way to prevent the effective even-channel chemical 
potential from becoming too large if the weakly-paired 
phase is to be established in the system. In this section, 
we propose a way to implement this requirement formally 
via generalized tunneling constraint. On the level of the 
BCS model used in Sec. HH this constraint leads to a 
stable weak-pairing phase in the even channel. In the 
following subsection A we describe a formal implementa- 
tion of the constraint in a system of BCS paired fermions. 
On the basis of this model, we propose, in subsection B 
a system in contact with a reservoir with which it can 
exchange particles so that, with the tunneling term in- 
cluded, a stable Pfaffian phase can be reached. 



A. A system with generalized tunneling constraint 

Cauchy determinant pairing describes p-wave pairing 
of j" and 1 particles, and in order to recover the spinlcss 
Pfaffian pairing, we need to "identify" j" with \.. Within 
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a functional formalism, that amounts to adding a term 
of the form 



X (r)[* t (r)-^(r)] 



(9) 



to the Langragian density via Grassmannian Lagrange 
multiplier x( r )- We will assume instead that we can al- 
ternatively express this via the constraint 

A(r)[M/t(r) - *J(r)][* t (r) - * ; (r)], (10) 

in terms of the bosonic multiplier A(r). By construc- 
tion this constraint affects only the odd channel. Within 
the mean-field approximation of a spatially constant mul- 
tiplier A(r) = A, one may identify A = Asas/%, i- c - 
the effect of the multiplier is the same as the tunnel- 
ing term in Sec. HH except for an overall decrease of 
the chemical potential, fx — > fx — A, which eventually 
yields a A-independent chemical potential in the even 
channel, fx e = fx, as mentioned above. Integration over 
the Lagrange multiplier projects to ^ ^ = 0, where 
^j, = [4'|(r) — vI>|(r)]/v / 2 is again the odd spin superpos- 
tion written in terms of the fermion fields vE , tT (r), i.e. it 
leaves us with no density in the odd channel. 

The BCS Hamiltonian including the constraint (fpQ|) 
has the same form as JT]) except that now ?k = ey — 
fx + A, as a consequence of the above-mentioned shift 
in the chemical potential. We can diagonalize it by a 
Bogoliubov transformation, 

a k = u-j-Ckt + uj.Cki + «T c -kt + v l c -ki- ( n ) 

The equation [ak, H) = Ea^ then defines the 
Bogoliubov-de Gennes equations and the Hamiltonian is 
transformed into the canonical form H = ^ k -E' Q k«k Q;k ~'~ 
Y^k EpkPkPk + E , where the eigenvalues ±E a k and 
±E/3k are given by 



E a -k = \/e k i + Agfc 2 , with e k ,i = e k - fx, 



Epk = J e k 2 + A 2 fc 2 , with e k , 2 = e k - fx + 2A, 

(12) 

respectively. The eigenvectors corresponding to E ak and 
Epu are, respectively, 



ak 



k y/E ak + e k ,i 
2\/ E ak k x iky 
A k 1 



(Ckt + c ki ) 
(cL kt + c f _ki) 



2\/^Qk \J E a k + 6k, 1 

k y/ £,3k + £ k .2 , . 
Pk = ~ r^— ~ — (Ckf-CkJ.) 



2\/Ef3k k x — iky 
A k 1 



2 y/ Epic \/^8 + e k,2 



(cltt-c^). (13) 



To find the stationary point for the action defined by the 
diagonalizcd BCS Hamiltonian at zero temperature, it is 



useful to continue A from the real axis to the complex 
plane C (see Ref. HH for details). The part Eq, through 
its dependence on A, 



■E 

k 



Agfc 2 



1 



E, 



ak 



£k.i 



E. 



pk + £k,2 



(14) 



determines the stationary point: dEg/dX = 0. Introduc- 
ing the notation fi — 2A = jl g C for the analytically con- 
tinued chemical potential in the odd channel, we continue 
as£/gk + ek,2 = —/t + o(fc 2 ), which gives us the station- 
ary point Re/i — > — oo or, equivalently, A — > +oo and we 
have strong coupling in the odd channel. For our choice 
fi > 0, we thus obtain weak coupling in the even chan- 
nel that has a Pfaffian description^ in Ckf + Ckj, (even) 
variables [see the expression for c*k (|13[) ]. Unlike the sit- 
uation in Sec. UH with the same assumption that Ao 
does not renormalize strongly with tunneling, the chem- 
ical potential of the even channel stays the same and the 
even channel can thus be weakly paired. 



B. A system with fluctuating density 

The previous discussion was on the simplified model of 
the 331 physics in terms of neutral fermions with an addi- 
tional constraint that leads to the Pfaffian, but seems ar- 
tificial and hard to implement in an experimental setting. 
Nevertheless it suggests a possible way to achieve the sta- 
ble Pfaffian phase. With respect to ordinary tunneling 
the generalized constraint can be modeled by strong tun- 
neling and an additional term in the effective description, 
AN, where N is the total number of particles and A tun- 
neling strength as before^ The chemical potential (of 
the whole system) depends on the tunneling and changes 
as fx — A. Then we have the following physical picture in 
mind: as we take A > this decrease of the chemical po- 
tential with tunneling will imply the decrease of the den- 
sity of the system. On the other hand, from the solution 
of the BCS system with the generalized constraint, we see 
that the effective chemical potential of the even channel 
stays the same - equal to fi, Eq. (fT2j) . This means also 
that the number of particles in the even channel stays the 
same, so the effects of the tunneling and the additional 
term cancels, but the polarization P increases with tun- 
neling. Thus we effectively maintain the same effective 
parameter fi with tunneling, its value will not increase, 
and we will be able to achieve the stable Pfaffian phase. 
To be more specific and quantitative about the role of 
fluctuating density in achieving the Pfaffian, we discuss 
in the remainder of this section the necessary dependence 
of density on tunneling. We will recover the demand for 
the decreasing density for the case of strong tunneling. 

Our BCS Hamiltonian, which is also the thermody- 
namic potential f2 at zero temperature, is specific for the 
fact that the independent thermodynamic variable, along 
with volume and temperature, is the chemical potential 
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/Lt A , given by /i A = \x — X. Therefore, as we change A (the 
parameter of the Hamiltonian) , we also induce a change 
in the chemical potential Sfi x = — X and this implies 



on 

~dX 



dE 
dX 



N, 



(15) 



for the particular system. The effective chemical poten- 
tial, on the other hand, is in this case 



eff 



/iP + (/i - 2A)(1 - P) = fx + (P - 1) 2A. (16) 



As before, this relation shows how the chemical potential 
renormalizcs as a parameter A of the BCS Hamiltonian 
is varied. We keep the volume constant and measure N 
with k F in the usual ansatz: 



N{k F ) = £ 



2 x 2tt 



(2tt)s 



dkk 



4/27T, (17) 



and Eq. (fT5"j) becomes 



N(k F ) = 



dE a (k Fl X)dk F dE (k F ,X) 



dk f 



dX 



dX 



(18) 



Differentiating and converting the sum in Eq into an in- 
tegral over k, we obtain the partial differential equation, 



k F 



A 2 k 3 



1 



1 



E a i 



dk k< 



e k F ,l 

A 2 



E, 



/3k F + ekj7,2 



S/Jk(£ ; /3k + £k,2) ' 



dk F 
~dX 

(19) 



from which we can extract the limiting case A ^> 
ti 2 k F /2m* when 



dN(kj 
dX 



= -cN{k F ), 



(20) 



with the constant factor c > 0, i.e. for large tunneling 
we should decrease the density of the system to stay at 
small [f and to stabilize the Pfaffian. We expect this 
limit to be pertinent to the experiments such as Ref. [25| 
where the compressible state starts to show signs of in- 
comprcssibility upon density imbalance. Enforcing the 
condition ([2TJ]) is then expected to lead to strengthening 
of the paired Pfaffian state, as summarized in Fig. [7] 
where the dashed line represents schematically a phase 
trajectory of the system that evolves under the general- 
ized tunneling constraint (in the large As^s limit). As 
Eq. (|2T)|) shows, in this case the density of the system 
needs to be decreased simultaneously with the increase 
of tunneling Ag^s . This can be achieved conventionally 
via the application of a gate voltage to the evaporated 
top/bottom gates as e.g. in Refs. UMM or by growing in 
situ a n + -GaAs layer that can serve as a gate4^ 

The main reason why we find it desirable to have an 
open system in the experiment is our inability to specify 
the parameters of the simplified, effective Read-Green 



model and pin-point the optimal density for which the 
Pfaffian state would be strong. With a better knowledge 
of the microscopic details of the system, one may as well 
fix a particular density (which we expect to be small) for 
a given, not too large, tunneling strength where the Pfaf- 
fian will be stable. However, an open system would allow 
a systematic study over a range of tunneling strengths in 
which the Pfaffian would show a characteristic strength- 
ening that would distinguish it from any other incom- 
pressible candidates. 




FIG. 7: Generalized tunneling contraint (dashed line) in the 
phase diagram of Fig. [4] Note that, as in Fig. 21 n denotes the 
effective chemical potential equal to the chemical potential of 
the even channel in the strong tunneling limit. 



Therefore, in principle, by changing the density of the 
system we can achieve a stable Pfaffian phase. We would 
like to compare the present discussion which is based on 
the simplified model of neutral fermions (mean field in 
nature with the simplifying assumption A Q = — A e = 
constant i.e. independent of tunneling) with exact di- 
agonalization in the LLL in the previous section. The 
arguments presented here call for an open system with 
adjustable density which demands also the adjustment 
of the magnetic field B to achieve the stable Pfaffian 
phase. Although doing so will preserve the filling factor, 
in general changing the total density may enhance the 
role of LL mixing and thus invalidate the LLL assump- 
tion in the exact diagonalizations. It will also lead to the 
renormalization of the parameters of the BCS effective 
model not taken into account in Sec. [TTJ Indeed, when 
lowering the total chemical potential via the generalized 
constraint, the density is also decreased. On the other 
hand the ratio between interaction strength (e 2 /Is) and 
cyclotron frequency (lo c ~ l/l 2 B ) is proportional to 1/ \fp 1 
from which we see that the LLL projection is invalidated 
if the density is significantly reduced. Therefore to reach 
and establish the Pfaffian phase it is likely that LL mix- 
ing has to be taken into account. This has been discussed 
in the recent literatur o 42 ' 43 as a way to stabilize the Pfaf- 
fian phase. Here we seek the Pfaffian in a two-component 
setting when a parameter of the system A is varied, which 
makes the inclusion of higher LLs harder. If we remain 
in the LLL, changing of the density amounts to simple 
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rcscaling of spectra (e 2 /lg — > c e 2 jls with c > 0), which 
cannot induce any signficant effect such as the change in 
the nature of (quasi)dcgcncracy of ground states on the 
torus. Even if the evidence for a Pfaffian phase is rather 
weak for the system sizes studied here, one may hope that 
the increase of these sizes will improve the case for such 
a state in the LLL as the odd channel may assume the 
role of the first excited Landau level before a complete 
polarization. 

V. CONCLUSIONS 

We investigated the possibility of creating the Moore- 
Read Pfaffian out of the paired two-component 331 state 
via tunneling. Exact diagonalization, performed under 
the constraint of LLL projection and fixed total number 
of particles, could not detect a stable Pfaffian phase, but 
a critical one between 331 and CFL phases. While the 
short-range interaction is likely to favor a direct transi- 
tion from the 331 to the CFL phase, long-range Coulomb 
interactions leave the possibility for a Pfaffian-like phase 
if the parameters of the system are tuned in a special 
way. Based on the connection between our numerical re- 
sults and the effective BCS Hamiltonian theory of paired 
states, we argue that one way to stabilize the Pfaffian 



state is to change the density (number of particles) of 
the system while increasing the tunneling. 

Our analysis was restricted to the Coulomb bilayer 
system and the tunneling term of the form —Asas^xj 
which is small in magnitude and generally difficult to con- 
trol. Similar considerations apply to the quantum well 
systems 2 ^ where f, J, stand for the two lowest electronic 
subbands and Asas acts like a Zeeman energy. In these 
systems, the analogue of the tunneling term used in our 
paper can be deployed to create asymmetric charge dis- 
tribution in the wide quantum welli^ The interplay of 
these two kinds of terms, tunneling and Zeeman, will be 
addressed in future work. 
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